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1. INTRODUCTION 



The geometric phase is now a standard topic in quantum mechanics [12] . 
To give an idea of what the geometric phase is, let us assume we have a 
parameter dependent Hamiltonian H(P(t)) where the parameters are time 
dependent [e.g., a time dependent magnetic field]. As a quantum system 
evolves under the Hamiltonian, it will pick up a total phase factor. We can 
decompose this phase into two parts. The first part is the dynamic phase 
which reflects how long the system has been evolving. The second part, the 
geometric phase part of the total phase, depends only on the path that P(t) 
takes in its evolution. That is, the geometric phase reflects the geometry of 
P's evolution in parameter space. An important point is that the geometric 
phase does not depend on how fast P(i) traverses its path, but depends only 
on the path itself. The article by Mukunda and Simon offers perhaps 
the most easily understood comprehensive introduction to the theory. 

Here we examine the dynamic and geometric phases in Hestenes' space- 
time algebra formalism This is an alternative mathematical treatment of 
the Dirac equation. So, the system we have in mind is a relativistic spin- 1/2 
particle evolving under a time dependent Hamiltonian. The time dependence 
of the Hamiltonian will come from a time varying electromagnetic field. 

The paper is organized as follows. In Section 2, we review the spacetime 
algebra formalism. The dynamic and geometric phase formulas are derived 
in Section 3. A discussion of the results follows in Section 4. 



2. SPACETIME ALGEBRA 

We review the spacetime algebra formalism that will be needed in the 
sequel. The reader is referred to for a more detailed account. We will let 
\l/ and $ denote the traditional [i.e., complex] wavef unctions, and ip and <p 
denote their spacetime algebra representations. 

Hestenes' spacetime algebra is the geometric [or Clifford] algebra of flat 
Minkowski spacetime. Before looking at this algebra, we will examine a 
simpler geometric algebra, the Pauli algebra. First, we will take as an axiom 
that the vector [or Clifford] product of a vector v with itself is given by 
vv := v ■ v, where v ■ v is the usual inner product. For our example we will 
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have three vectors 01, cr 2 and 03. Let v = o~\ + cr 2 so that 
vv = (a i + er 2 )(<Xi + cr 2 ) 

= CTi<Ti + (Ti^ + (T 2 (Ti + <7 2 cr 2 (1) 

and 

V • V = (<7l + CT 2 ) ■ (o"l + (7 2 ) 

= (7l • (7l + 2(Ti • (T 2 + (7 2 • CT 2 . (2) 

Since vv = v ■ v, and using 0\0\ — G\ ■ o~\ and <7 2 <7 2 = ct 2 • cr 2 , we see from 
© and © that 

0"! • cr 2 = ^(o"io- 2 + a 2 cr 1 ). (3) 

Equation Q allows us to define the inner product in terms of the vector 
product 

a i -a j = -(a^j + a^i). (4) 

Because a^-aj is a scalar, we can form a metric tensor by defining rjij := <Ji-aj, 
where i,j = 1,2,3. Then (@J becomes 

Oi ■ <Tj = Tjij. (5) 

Now notice that 

(Tj = ~{(JiOj + OjOi) + ~{(Ji<Jj - OjGi). (6) 



a 



The first part of (jH)) is the inner product of o~i and a r The second part of 
(JUJ) we will call the outer product of <jj and o~j, and denote this by o~i A o~j := 
l/2(crj<7j — OjCri). So (0) can be rewritten as 

OiOj = (Ti ■ o~j + (Tj A Oj. (7) 

What happens if cij = Aoj, where A is a scalar? Then, 

(TiAo-j = ^ApiOj - o-j(Ji) (8) 



A 

2 
0. 
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So, if <7j and Oj are not a multiple of each other [i.e., if they are not collinear], 
our vector product will have a nonzero outer product. 

The geometry of geometric algebra is seen by how we interpret all of this. 
As usual, we consider the inner product cr, • Oj as being the projection of <7j 
onto Oy But what of cr, A op. This we let represent the area between the 
vectors cr, and Oj in a plane that contains both vectors. 

Let us now make some refinements. First, instead of representing a [real] 
scalar as simply A, we will denote it by Al. Since scalars will commute with 
vectors [e.g., \o\ = oi\], we will have 

l<7i = ad. (9) 

Now, what if our vectors are orthonormal so that 

<7j • (Tj = 5ij7 (10) 

Then, if % ^ j, 

OiOj = di ■ a j + (Tj A a 'j 

= (Tj A Oj 

= —Oj A Oi 

= -OjOi. (11) 

So, if Oi and <7j are different orthonormal vectors, the outer product <7j A (Tj 
is not simply an area. Rather, it is a directed area since OiOj = —OjOi. So 
it depends on whether we go "clockwise" or "anti-clockwise" over the area. 
This is a generalization of the concept of a vector's direction. To formalize 
this new concept, call the elements o^j bivectors. Since i,j = 1,2,3, there 
are three independent bivectors given by <7i<7 2 , <7i<7 3 and o 2 o 3 [note that, for 
example, &i0 2 = —o 2 0\ and o 3 o 3 = o 3 ■ o 3 ]. A further generalization is to 
have a directed volume element, or a trivector, X := o\0 2 o 3 . The trivector 
satisfies 

J 2 = (<7i<7 2 <7 3 )((7i<72<73) 
= -<7i<7i<7 2 <7 2 <73<73 

= -I- (12) 

So, by starting with three vectors <7i, 02 and o" 3 and a vector product, 
we are naturally lead to new elements 1, X and the bivectors. That is, 
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given three vectors, we can imagine, in addition to the vectors themselves, 
a scalar component, a directed volume and directed areas. So our basis 
for this construction is given by {1, Oj, o,,Ofc, 01O2O3}, where i = 1,2,3 and 
1 < j < k < 3. Notice that 



0io 2 3 0i 

0j0102O 3 
0203 



(13) 



and that Xa 2 = —0103 and X03 = 01 02. So we can also use {1, o^Xo^X}, 
z = 1, 2, 3, as our basis. We prefer to use this latter basis. 

It is possible to represent our vectors as matrices. For our example, the 
familiar Pauli spin matrices 



01 



" 


1 " 




" -i ' 




" 1 


1 





02 = 


i 


£3 = 


-1 



(14) 



where i = \f— T, are one possible representation. [In a matrix representation 
we let 1 be the identity matrix.] While it is nice to have an explicit matrix 
representation for our vectors, we must stress that they are to be treated 
as vectors. That means that under a spatial rotation, the basis vectors will 
change. For example, if we rotate in the 0102-plane by an angle (p, a vector 
will transform as 



w 



w = UwU 



where U = exp[— Xo 3 (^/2]. So, if (p — ir/2, we expect that a\ 
and 03 — > 03. This is easily checked since, e.g., 



02, 0"2 



(15) 
-01 



exp 



-X03 — 
4. 



02 exp 



7T 



X03 — 

4 



exp 



— Xo s — 



cos 



7T 

2 

X0302 



02 

Xo 3 sin 



02 



= -01, (16) 

where we used the facts that Xoj = 0, ; X in the first line, and that (Xoj) 2 = — 1 
in the second. Equation (J15)) also holds for a general rotation using the Euler 
angles (ip, 8, x)- In this case U has the form [3] 



U 



exp 



exp 



-X05- 



exp 



-X03- 



(17) 



5 



So far we have only considered scalars, vectors, etc. We can also have 
more general objects, called multivectors, which are of the form 

C = c s + (c> 1 + c^ 2 + c^ 3 )+X(cJa 1 + 4a 2 + 4a 3 ) + c'J, (18) 

where C has scalar [c s ], vector [c"], bivector [c b ] and trivector [c*] parts. Let 
us introduce two operations on C that will prove useful in the sequel. The 
first is to let (C) s denote the scalar part of (fTHjl. so (C) s = c s . For higher 
parts of C, we will let (C)i be the part of C that can be minimally expressed 
using i vectors. So, (C)i = c" x a\ + c 2 cr 2 + C3V3, (C) 2 = 1(c\a\ + c 2 cr 2 + C303) 
and (C)s = c l 1. The second operation is the reversion of C, written as C. 
By this we mean the reversal of all the vector products in C. So, for example, 
<j\G 2 = o~i<J\. In general, for two multivectors A and B, we have AB = BA. 
For (fTHj) notice that X = —X and XcTj = — a(L = —Xoi [because XcTj = a(I\ ■ 
Then 

C = c s + (c>i + c v 2 a 2 + c v 3 a 3 ) - l{c\a x + c b 2 a 2 + c b 3 a 3 ) - c% (19) 

We can use the reversion operation to rewrite p5jl . From (fTTjl we see that 
U^ 1 = U. So (fT5|) can be written as 

c -> d = UCU, (20) 

where we now allow the rotation to act on a general multivector. 

Now we will review Hestenes' spacetime algebra. The spacetime algebra 
is generated by four vectors {7^}, fj, — 0,1, 2, 3, that satisfy the Dirac algebra 

7^-7^ = 9p, (21) 
= diag(H ). 

The metric tensor g^ v in (|2*Tj) is the usual one from special relativity. The 
7o vector is the time direction and the j n , n = 1,2, 3, vectors are the spatial 
directions, in an observer's frame of reference. We can also raise the indices 
on the 7^'s by defining the vectors r y u as those that satisfy the relation 7^-7^ = 
5". [Notice that 7 = 70 and 7™ = — 7„ satisfy this.] One possible matrix 
representation of the 7 M 's are the standard Dirac matrices. 
The vectors 7 M result in the basis 

{1,7^, (cr n ,Xcr n ),X7 A1 ,X} , (22) 
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where n = 1, 2, 3 and a n : = 7„7 . Also, X := 7o7i7273 = ^l^a, typ = -7 M X 
and now X = X. The Pauli algebra we examined above is a sub-algebra of 
the Dirac algebra. So the 1, a n , Ta n and X in (J22j) satisfy the geometric 
algebra we looked at using the <7 n 's as our vectors. We can consider the 
Pauli algebra as the algebra of space and the Dirac algebra as the algebra 
of spacetime. Note that in the Pauli algebra, the <r n 's are vectors while the 
Xcr n 's are bivectors. In the Dirac algebra, both of these are bivectors. 

In the spacetime algebra formalism, a proper, orthochronous Lorentz 
transformation [i.e., a spacetime rotation that does not reverse time or space] 
is represented by R. Since every Lorentz transformation can be decomposed 
into a pure boost and a pure spatial rotation, R will have the form R = LU, 
where U is the pure spatial rotation given in (JT7J), and the pure boost, L, 
has the general form 

(bxax + b 2 o 2 + 630-3) 
L = exp 

The spacetime generalization of is given by 

w -> w' = RwR- 1 . (24) 

We have already seen that U^ 1 = U which still holds in the Dirac algebra 
[because a n = —<r n , X = X and Ia n = cr n X, so (Xcr n ) = — Xcr n ]. Because 
o n = —cr n , we also have that L _1 = L, so R = UL = R" 1 . As we did in (|20|). 
we can write (|24|) as 

C -> d = RCR, (25) 

where C is a general spacetime multivector. 

It can be shown [Sj that a Dirac spinor, in the spacetime algebra, has the 
general form 

V = (pe^) 1/2 i?, (26) 

where p(x) and f3{x) are scalars and R(x) is a [proper, orthochronous] Lorentz 
transformation. [Deriving ()26|) is easy though lengthy, so we will simply 
accept this form for ip without proof.] Note that p, j3 and R depend on 
the spacetime point x. That is, they are locally defined parameters in the 
wavefunction ip. In particular, R is a local Lorentz transformation- not a 



(23) 
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global one. This observation leads to the local observables theory favored by 
Hestenes jjj. Now, the p term in is interpreted as a proper probability 
density. A satisfactory interpretation of the j3 parameter is still lacking in 
Hestenes' theory, however. So Hestenes' interpretation of the Dirac equation, 
while attractive, is not yet complete. Using the form for R above, can 
be rewritten as 

if, = [ pe ^) 1/2 LU (27) 
= (pe^LUoe- 1 ^ 2 , (28) 

where, from (fT7f). 

Multiplying \I/ in the standard Dirac theory by a phase factor exp[ia] 
corresponds, in the Hestenes-Dirac theory, to 

ip -> iff' = ipe Iasa . (30) 

Equation (|50|) simply expresses a correspondence between two different math- 
ematical formalism, the standard Dirac theory and Hestenes' theory. It is 
useful to note two things though. First, Ta 3 , in Hestenes' formalism, corre- 
sponds to the complex number i in the standard Dirac theory. Also, given 
(J28|) . multiplication by a phase factor in the standard Dirac theory corre- 
sponds to a change in the Euler angle \ in the spacetime algebra formalism. 
This allows us to identify — x/2 as the phase in Hestenes' theory. [Note the 
minus sign here. In (|30|) we added a new phase a to the original phase. The 
minus sign follows from (|28|).] Another lengthy, but easy, derivation shows 
that the Hestenes form of the Dirac equation for a spin- 1/2 particle of charge 
e and mass m is given by [with H = 1 and c = 1] 

OipXa 3 — eAip = rmp'jQ, (31) 

where A := is the electromagnetic field vector and □ := 7^9^. [The 

Einstein summation convention is used here, so J^A^ := 7°v4o + 7 1 ^4i + 
7 2 v4 2 + 7 3 A 3 .] 

In the standard Dirac theory, the probability density g in an observer's 
frame of reference is associated with the time component of the Dirac current 
ty^ty, where is the Hermitian adjoint of the complex wavefunction In 
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order to translate this into Hestenes' formalism, we need to find the spacetime 
algebra equivalent of where $ is some other complex wavefunction. The 
spacetime algebra representation of is given by Q 

tft ^V f := 7o^7o- (32) 

Now, = r + ic where r and c are real scalars. The real part of 

ty^<&, denoted is simply the real scalar r. This corresponds to the 

scalar part of (pjoi^lo, the spacetime algebra representation of Then 
sft(\j/t$) <_> (07qV'7o}s- Now recall that i in the Dirac theory corresponds 
to Zo" 3 in Hestenes' theory. Then the imaginary part of W<f>, written as 
= — goes to — (0Xo"37o , ?/ ; 7o)s in the spacetime algebra for- 

malism. So the spacetime algebra representation of the complex probability 
amplitude density of standard quantum theory is given by [T] 

(07oV>7o) s -Xo- 3 (0X0-370^70)5- (33) 

Equation ()33|) is equal to the probability density g when $ = 

Notice the ubiquitous presence of 70 in (|3*2|) and (JBlIj) . Since 70 is the time 
vector in an observer's frame, we say that the Hermitian adjoint W and the 
probability density g are "frame dependent". That is, they depend on the 
direction of 70, and, hence, on the observer's frame of reference. 

Results similar to (|3~2|l and (|3*3*|l hold for the Dirac adjoint \l/ := 
where 70 is the standard Dirac matrix. The spacetime algebra representation 
of ^ is P 

vF «-> <j>. (34) 

Also, 

<-> {^) a -Xa^Xo^) s . (35) 

Since f)34jl and (J3*5|) do not contain a 70 term, we say they are frame inde- 
pendent. 

The ip in (|26|) can also be rewritten in a form like that in (|18|) . That is, 
ip can be expressed as 

ij) = c s + + <%a 2 + c v 3 a 3 )+l{c b 1 a 1 + c b 2 a 2 + c b 3 a 3 ) + c* J, (36) 
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where a n = 7 n 7 in the Dirac, versus Pauli, algebra. Now, all the c's in ()36|) 
are real, so it would appear that we have managed to remove the complex 
structure of the Dirac theory by using the spacetime algebra. This is incor- 
rect, as (}3"3"j) shows. The "1" and "Xo" 3 " of Hestenes' theory correspond to 
the real and imaginary parts, respectively, of the standard Dirac theory. So, 
even though we can express if) completely in terms of real numbers, as in 
(I36j) . the complex nature of the Dirac theory still remains. As explained in 
such a complex structure is inherent to any correct formulation of Dirac's 
equation. 

Lastly, we need to define a few terms that will be used latter. [To rig- 
orously justify these definitions will take us too far afield. The reader can 
consult jSj, and the references therein, to find their physical justifications. 
For now, they can simply be taken as definitions.] The proper velocity vec- 
tor v = v is defined by 

v := Rj R (37) 
= #o7o£o, (38) 

where R is from the wavefunction if) and Rq := LUq. That is, v is given by the 
spacetime rotation of the time axis 70 determined locally by the wavefunction 
if). We have that vv = 1, so the velocity vector is normalized. From (|33|) we 
can see that, when $ = 

q = pv° } (39) 

where p is also from the wavefunction if). We define the angular velocities as 

n„ := 2{d tx R)R (40) 
^ := 2{d fl R )R . (41) 

Because RR = 1, we have that (d fl R)R = —Rld^R). Notice that this implies 
that = — O^. Similarly, we can show that uj^ = —uj^. Now, the spin 
polarization vector s = s^j^ is given by 

s := -R l3 R (42) 
= |Ro7 3j Ro. (43) 
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Finally, the spin angular momentum bivector S = Xsv is defined by 

S := ~RIa 3 R (44) 

= ^R Ia 3 R . (45) 

Equation (jUj) gives us d^S = l/2(Q fl S — SQ^). Similarly, from we have 
d^S = l/2(o; M 5 — Su^). From (JHHJ), (|4l?jl and (|43|) we see that i?o determines 
the orientations of v , s and S 1 . By orientation we mean the directions of the 
vectors v and s, and the "tilt" and direction of the directed area S, in the 
four dimensional spacetime. 



3. PHASE FORMULAS 

Here we derive the spacetime algebra formulas for the dynamic and ge- 
ometric phases. Initially we simply translate the standard dynamic phase 
formula into Hestenes' formalism. But, much of Hestenes' theory deals with 
local observables jjj, so we will define a local dynamic phase. This local 
definition is useful for working in the spacetime algebra. Next we derive the 
geometric phase formula. We then briefly examine the relationships between 
the spacetime algebra and standard formulas. Finally, we will find that we 
are able to redefine the phases, allowing us to use much simpler formulas in 
Hestenes' theory. 

First, notice that (jHljl is invariant under the transformation 

ip — > ip ex P P^s^L (46) 

where a is a real constant. In the Dirac theory, this corresponds to adding a 
constant phase to In Hestenes' theory, we are adding a constant angle to 
the Euler angle x i n (EEJ)- Let us now define the projection operator n by 

n(^) := ^tp : ip = ipe Xa ' Aa } for all real constants a\ (47) 

and let V be the space of all such projections. This amounts to projecting 
the multivectors tp in the Hilbert space TC onto the representative ray n(^) in 
V, where Tl(ip) differs from ip only by a constant spatial rotation exp[Xo"3a]. 
[Equation (}4Tj) corresponds to the projection of \& onto the ray IL(\I/) , where 
n(\I/) differs from \1/ only by a phase factor exp[ia].] If VKO evolves along the 
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curve C in TC, then II(^>(£)) will evolve along the curve C in V. Notice that 
we have parameterized ?/>'s evolution in Ji by £. The geometric phase should 
only depend on the curve C in V, the path of n(^) in V. In particular, the 
geometric phase needs to be independent of the rate at which II("0) traverses 
C. Therefore, a reparameterization of £ can not affect the geometric phase. 
The [complex] global dynamic phase, A@, is given by 

A G = J d£d 3 x S(tfttf), (48) 

where the overdot represents differentiation with respect to £. As we have 
seen earlier, ^(SIZ+SI/) <->• — (ipZa^Q^o) s . So the spacetime algebra formula 
for the global dynamic phase, 5g, is 

5 G = - J d£d\ (^Ta 3 7o^7o) s - (49) 

From (|49|) let us define the local dynamic phase by 

q h : = -(^Jo- 3 7o^7o) s . (50) 

We use g, rather than p, in (J5Uj) because the Hermitian adjoint used in (J4*9~j) 
is frame dependent, so we expect the same for our probability distribution. 
Also, since the Hermitian adjoint singles out a preferred time direction, we 
will let £ be the observer's time, given by t. Using (|27)|) in (f57Ij) results in, 
after some tedious algebra, 

q 5 L = pP(- R^R^s - p{RT^Rlo)s 

= p$(sj ) s - p(fi T,s7o) s 
= ps°p -p(Qq Svj ) s 

= g^P-p{n S(v° + v)) s , (51) 

where vy := v° + v, v := v n a n , and Sv = Xs because vv = 1. The Q S term 
in (|BT|) can have scalar, bivector and pseudoscalar [i.e., X] parts, denoted by 
Qo ■ S, (QoS)2 and (floS)^, respectively. Under the reversion operation, the 
scalar and pseudoscalar parts are even, while the bivector part is odd. Hence, 

(n s) 2 = ±(n s-(^)) 

= ^(n Q s-sn ). (52) 
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From (J44j) . we see that 



Using this in (|5T|) gives us 



-{Rla 3 R + Rla 3 R) 

(n s) 2 . 



(53) 



s° ■ 

= Q ~n$ - Q^o ■ S -p(Sv) s . 



(54) 



Now, s ■ v = 0, so 



s ■ v 



(55) 



where 370 = s° + s. Thus, (p>4"j) becomes 



-o n • s - — • 



From pity to (jHHjl . it follows that, after letting £ = t, 



>G 



dtd 3 x gSi 



dtd x q \ —Qo ■ S • 



(56) 



(57) 



We now find a spacetime algebra geometric phase formula. First, let us 
define a new wavef unction ip := ip exp[— I<j 3 5l] that differs from ip only in 
having the local dynamic phase removed. Now, 



4 



+ ipe 



-Ia 3 5 L 



(58) 



Using (J3l?j) for the inner product of ip' with ip' , and noting that (^Zo^oV'To).' 
0, we have from (1581) 



(^7o^7o) s - Jo- 3 (^Jo- 3 7 ^7 ) s 

-Jo- 3 5l(V ; 7o^7o)s 

(?/>7o^7o) s + Icr 3 g5 L - lcr 3 g5 L 

(hoho)s- (59) 
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Since (4>' lo 1 ^' lo) s — (V'7o'07o) s, it follows from (JoT?|) that 

(^'Ta 3 7o^7o)s = 0. (60) 
With ip' differing from ip only by the x factor of (|28|). let us write 



(pe X/3 ) 1/2 Lf/ e- Iff «' /2 . 



(61) 



Using (|5T|) in (|5Uj) . and doing some lengthy algebra, gives us, similar to our 
derivation of (J53I) . 



(^'Xa 3 7o^'7o) ; 



!<(-i2oX<7 3 x + to + 2i? )Xa 3 7o J Ro7o), 
£> , s° • 

77X - £^t/? + ^0 • S + p(Sv) s 
2 v u 

g ■' s ° h , c , v c 

77X - Q-kP + Quj ■ S + g— ■ S, 

2 v u ir 



(62) 



where we used the fact that S = (co , S') 2 [see fB3j) ]. Then, from (|50|) and 
we have 



2 X 



-^0-5' 



(63) 



Notice that all of the terms on the righthand-side of (J63|) depend only 
on the path C in V. That is, they are determined by R®, not by R. [This 
does not hold for 6l because of the presence of fin in (|5B)).] Also, both 
sides of (J53j) are linear in the time derivative. Hence, \ is independent of a 
reparameterization of t. These two observations allow us to define the local 
spacetime algebra geometric phase by jl := — x /2 + x(0)/2. [We use a 
minus sign here because — x /2 = — x/2 — 5l- The x(0)/2 factor is because 
7l(0) must vanish. Then 7l + 5l = — x/2 + x(0)/2, where we have previously 
identified — x/2 as the total phase.] Thus, 



7l 



-fin • S ~ 



UJ ■ S + 



5 



s •" 

n0 , 



(64) 
(65) 
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where <5l(0) = and 7l(0) = 0. To find the global spacetime algebra geo- 
metric phase, 7g, we need to perform the integration 

JG •= J dtd 3 X Q'Jl 

= J dtd 3 x g j^o • S + ■ S - I , (66) 

similar to ()57|) above. 

Equation (JMj) is a straightforward translation of the standard dynamic 
phase formula into the spacetime algebra. Thus, the value of Sq will equal 
that of Aq- We may ask if the same holds true for jg an d the standard 
geometric phase Tg? This will not generally be the case. The — x/2 factor in 
()28|). which we identified as the total phase, may contain dynamics beyond 
the local representations of Aq and T^. As an example, let us consider an 
adiabatic evolution of a wavefunction that starts in an energy eigenstate. Let 
the m-th energy eigenfunction be given by 

K = [pe^f 2 R m0 e-^'\ (67) 

The adiabatic theorem states that our wavefunction ip m can differ from m 
only by a phase factor. Thus, using the complex global phases A m G and T m G 

^ m = m e I<T3(A - G+r - G) . (68) 

[Strictly speaking, (J68j) will contain some 0(e) terms where e is an infinites- 
imal number that reflects the degree of adiabaticity, and the 0(e) terms are 
due to departures from strict adiabaticity. We will ignore these terms in 
(ff)5|) .] In this situation, T mG is given by [12] 

r mG = -J d 3 x% ($J n $ m ) . (69) 
A derivation similar to that in Section 3 for the dynamic phase shows that 

= ImG + J d 3 X^. (70) 

This result does not invalidate calling 7^ the spacetime algebra geometric 
phase. It still represents that part of the total phase — x/2 due only to the 
geometry of the wavef unction's evolution. 
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Now, (|26|) can be rewritten as 



$ = (pe J/3 ) 1/2 J R e lCT3(7i+<5i - x(0)/2) 



(71) 



Thus, we can redefine the phases as 



5 L 




(72) 
(73) 



where 5l(0) = and 7l(0) = 0. Equations (J?2J) and (J73*|) are still physically 
meaningful definitions. They capture all the essential properties of the dy- 
namic and geometric phases. The simplicity of (|72)) and (J75j) . versus (|54*Jl and 
respectively, suggests that they are more useful definitions for Hestenes' 
formalism. It is easy to show, for example, that ()72j) and ()73|) hold exactly 
in the non-relativistic limit. That is, the second terms on the right-hand 
sides of ()64j) and (J55j) vanish in the non-relativistic limit because v — > 0. 
Thus, we can use exactly the same phase formulas in the relativistic and 
non-relativistic cases. Also, it is easy to show the (|64j) and (|65j) reduce to 
(J72J) and (|73*|). respectively, when we consider the adiabatic evolution of an 
energy eigenstate. Overall, (|T2"j) and (J73*|) seem to be more reasonable defini- 
tions of the phases, in the spacetime algebra, than ()64|) and (|65|). Finally, we 
note that Hestenes' had previously proposed using Q ■ S as the total local 
phase formula ^U]- As we see, except for the sign, this is only the dynamic 
part of the phase. We also need to take into account the geometric part of 
the total local phase. 

The formulas (|T2"j) and (|73*jl can be written as 



Unlike (|32j) and (J33|) . these no longer have 70 present. Thus, they are frame 
independent quantities. Now, the wavefunction if) determines a set of stream- 
lines via the proper velocity vector v. That is, if we imagine a particle as 
actually starting at some initial spacetime point xq, the spacetime vector v 
will determine its future positions, given by the streamlines. Along a given 
streamline there is a proper time r. Because (|74*j) and ({75} are independent 



-(KZa 3 R) s 
(R 2cr 3 R ) s 



(74) 
(75) 
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of a particular reference frame, we can allow the derivative to be with respect 
to r. Hence, (J74*|) and (|75j) can also be used for the proper phase formulas 
[remembering that, in this case, we are using a proper time derivative]. 



4. DISCUSSION 

Let us first review a few facts about (3 that will be useful in the fol- 
lowing discussion [see the article by Gull, Lasenby and Doran [S] for more 
background]. When there is no electromagnetic field, i.e., A = in ([31)1 . 
the Dirac equation admits plane wave solutions. For the electron solutions 
(3 = 0, while for the positron solutions (3 = tt. However, when A ^ 0, a 
general wavefunction can have other values of (3, as demonstrated by the 
[non-relativistic] solutions for the hydrogen atom. Additionally, in their nu- 
merical simulations of tunnelling times, Gull et al. show that p and j3 are not 
necessarily constant at a given position for all time. So the (3 term in 
does not necessarily vanish. 

Now, it can be shown that 

7 5 # <-> (76) 
pcos((3) (77) 
^75* <-> -psin(/?), (78) 

where 75 is the standard Dirac matrix. It follows that 

<-> pe X/3 . (79) 

So the e J/3 term can be thought of as a local chiral transformation |2J, an 
observation made previously Let us take this literally and think of (3 as 
the chiral angle in the Dirac wavefunction. 

Equation ()64|) for the dynamic phase is a straightforward translation of 
the standard phase formula ()48|) into the geometric algebra. We can rewrite 
as 



A G = J d 3 x S(^) 



- / d 3 x ^ ] i^m. (80) 
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In the standard Dirac theory, (|8U|) is interpreted as the negative of the ex- 
pected value of the energy operator. The same physical interpretation is then 
given to locally However, interpreting (JMj) using only the spacetime 

algebra formalism is difficult because of the presence of j3. 

In contrast, the dynamic phase formula in (|72j) is easily interpreted in the 
spacetime algebra. It is the negative of the component of Qq in the spacetime 
plane S. To translate (|72j) back into the standard Dirac theory, we use (I74j) 
and ()35|) . If R is the spacetime representation of R, in the standard theory 
d?U) is given by 

A L = Qf(RR). (81) 

Because RR = 1, we have from (J79|) that RR =1. It follows that RR = 
— (RR)^. Hence, RR is purely imaginary. So (JHTj) becomes 

A L = -iRR. (82) 
The question now is, how to find R? Since 

<-> lip (83) 
we can write ()26|) in the standard theory as 

m = ^e^ /2 R. (84) 
Provided p ^ 0, we have that 

R = ^. (85) 

VP 

So R is given by a local chiral transformation of Notice that (j84j) implies 
that R is invariant under a local chiral transformation of 

For both sets of phase formulas, we see that they are easily interpreted 
in one formalism but not the other. Both geometric phase definitions are 
invariant under a local gauge transformation. [This likely accounts for the 
difference in 7 G and T G , see (fTUj) .] The phases in (|T2"j) and (fT3~j) are also 
invariant under a local chiral transformation. These may prove to be more 
useful in the electroweak theory. 

It is difficult to decide theoretically which set of formulas is the correct 
one. It may also be difficult to experimentally verify which set is correct. 
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This is because (|64|) and (|B5j) reduce to (|72|) and (J75|) for adiabatic evolu- 
tions of energy eigenstates and, in the non-relativistic limit. 
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